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CN ■ Abstract 

\ We study the global existence of weak solutions to a multi-dimensional simplified Ericksen-Leslie 
pL| I system for compressible flows of nematic liquid crystals with large initial energy in a bounded do- 
CN ■ main Vt C M^, where = 2 or 3. By exploiting a maximum principle, Nirenberg's interpolation 
inequality and a smallness condition imposed on the A^-th component of initial direction field do 
to overcome the difficulties induced by the supercritical nonlinearity |Vdpd in the equations of 
angular momentum, and then adapting a modified three-dimensional approximation scheme and 
the weak convergence arguments for the compressible Navier-Stokes equations, we establish the 
global existence of weak solutions to the initial-boundary problem with large initial energy and 
without any smallness condition on the initial density and velocity. 
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■ 1. Introduction 

(N 

. We study the global existence of weak solutions to the following mult i- dimensional simplified 

version of the Ericksen-Leslie model in a bounded domain Vt C which describes the motion 
of a compressible flow of nematic liquid crystals: 



O 



+ div(pv) = 0, (1.1) 



^ ■ dt{p^) + div(pv (g) v) + VP = pAv + (/i + A)Vdivv 

-z/div(^Vd0 Vd- i|Vd|2lj, (1.2) 

9td + V- Vd = e(Ad+ iVdpd), (1.3) 

with initial conditions: 

p(x,0) = po(x), d(x,0) = do(x), (pv)(x,0) = mo(x) in fi, (1.4) 

and boundary conditions: 

n ■ Vd(x, t) = 0, v(x, t) = 0, xeOn, t> 0, (1.5) 
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where n denotes the outer normal vector of f2. The unknown function p is the density of the 
nematic hquid crystals, v the velocity and P{p) the pressure determined through the equations 
of state, d represents the macroscopic average of the nematic liquid crystal orientation field. 
The constants /i. A, i/, and 6 denote the shear viscosity, the bulk viscosity, the competition 
between kinetic and potential energies, and the microscopic elastic relation time for the molecular 
orientation field, respectively, they satisfy the physical conditions: 

/i>0, A + /i>0, z/ > 0, 9>0. 

I denotes the N x N identity matrix. The term Vd © Vd denotes the N x N matrix whose 
(i, j)-th entry is given by d^^d ■ d^jd, for I < i,j < N, i.e., Vd Vd = (Vd)'^Vd, where (Vd)'"' 
denotes the transpose of the N x N matrix Vd. 



In 1989, Lin [17| first derived a simplified Ericksen-Leslie system modeling liquid crystal flows 



when the fluid is incompressible and viscous. Subsequently, Lin and Liu [19|, |20| established some 
analysis results on the simplified Ericksen-Leslie system, such as the existence of weak and strong 
solutions and the partial regularity of suitable solutions, under the assumption that the liquid 
crystal director field is of varying length by Leslie's terminology, or variable degree of orientation 
by Ericksen's terminology. 

Since the supercritical nonlinearity |Vdpd causes significant mathematical difficulties, Lin in 
171] introduced a Ginzburg-Landau approximation of the simplified Ericksen-Leslie system, i.e., 
|Vdpd in fll.3p is replaced by the Ginzburg-Landau penalty function (1 — |dp)/e or by a more 
general penalty function. Consequently, by establishing some estimates to deal with the direction 
field and its coupling/interaction with the fluid variables, a number of results on the Navier-Stokes 
equations can be successfully generalized to such Ginzburg-Landau approximation model. For 



examples, when p is a constant, i.e., the homogeneous incompressible case, Lin and Liu [19 
proved the global existence of weak solutions in 2D and 3D. In particular, they also obtained 
the existence and uniqueness of global classical solutions either in 2D or in 3D for large fluid 
viscosity /i. In addition, the existence of weak solutions to the density-dependent incompressible 
flow of liquid crystals was proved in H, 23|. Recently, Wang and Yu 3^, and Liu and Qin 



2J] independently established the global existence of weak solutions to the three-dimensional 
compressible flow of liquid crystals with the Ginzburg-Landau penalty function. 

In the past a few years, progress has also been made on the analysis of the model fll.ip - 
(11.31) by overcoming the difficulty induced by the supercritical nonlinearity |Vdpd. For the 
incompressible case, the existence of large weak solutions in 2D was established in [isj and 



13( 1 for a bounded domain and the whole space respectively, and the local existence of large 
strongsolutions and global existence of small strong solutions in three dimensions were proved 
in pj, Ig, [l5|, l2l|, |29j. For the 3D compressible case, the existence of strong solutions have been 



investigated extensively. For examples, the local existence of strong solutions and a blow-up 
criterion were obtained in js], [sj , while the existence and uniqueness of global strong solutions to 
the Cauchy problem in critical Besov spaces were proved in [f\ provided that the initial data are 
close to an equilibrium state, and the global existence of classical solutions to the Cauchy problem 



was shown in ^J] with smooth initial data that has small energy but possibly large oscillations 
with possible vacuum and constant state as far-fleld condition. Recently progress has also been 
made on the existence of weak solutions to multi-dimensional problem fll.ll) - (ll.3p . For examples. 



Jiang et al [10| established the existence of global weak solutions to the two-dimensional problem 
in a bounded domain under a restriction imposed on the initial energy including the case of 
small initial energy. Moreover they also obtained the existence of global large weak solutions 
to the two-dimensional Cauchy problem, provided that the second component of initial data of 



2 



the direction field satisfies some geometric angle condition. At the same time, Wu and Tan 31 
established the existence of global weak solutions to the Cauchy problem fll.ip -( lTT^ by using 



Suen and Hoff's method [28|, if the initial energy around equilibrium state is sufficiently small, 
the coefficients fi and A satisfy < X + fx < (3 + \/21)fj./6, and the initial data (vq, do) satisfies 
||vo||lp(r3) + ||do||LP(R3) < oo with p> 6. 

To our best knowledge, however, there are no results available on weak solutions of the 
multi-dimensional problem f ll.ip -( lL3l) with large initial data in a bounded domain, due to the 



[])-(0 

difficulties induced by the compressibility and the supercritical nonlinearity. It seems that the 
only global existence of large weak solutions to fll.ip -( !L3|) was shown in the ID case in jij. 
On the other hand, there exists a global weak solution to the mult i- dimensional compressible 
Navier-Stokes equations with large initial data (i.e., the initial energy can be arbitrarily large). 
A question naturally arises whether one can establish a global existence result for the problem 
f ll.ip - fll.Sp without any smallness restriction imposed on the initial density and velocity. In the 
current paper, we give a positive answer to this question in the two-dimensional case under a 
restriction on the last component of initial direction field do, while in the three-dimensional case, 
a somewhat weaker existence result is obtained. 

Before stating our main result, we explain the notations and conventions used throughout 
this paper. In this paper we focus our study on the case of isentropic fiows as in j3oj and assume 
that 



For the sake of simplicity, we define 



with A > 0, 



/:=/T:=(0,r), QT = nxI, 



J^{t) := J^(p,v,d) := I (/i|Vv|2 + (A + /i)|divv|2 + 0(|Ad+|Vd|2dn)dx, 

Jo. 



;i.6) 



and 



^/ N ^/ ,N /■ Imp A ^ z/^|VdP\ ^ ^ ^ 

8[t) := £{p,m,d) := j ( 2^^^^^'"^°^ ^ 2 / ^ m = pv, (1.7) 

where l{p>o} denotes the characteristic function. We use the bold fonts to denote the product 
spaces, for examples, 

:= {LP{n)f, H^(fi) := {H^{n)f = {Wl;'\n)f, H^(fi) := {w''^\n)f; 

and the Sobolev space with weak topology is defined as 

C°(/, Ke.^.m := {f : / -> ^'{^) / f • ^ g G L^(fi)| . 

In what follows, the letter Cq will denote a generic positive constant which may depend on the 
dimension of space N, and the letter C(. . .) will denote a generic positive constant depending on 
its variables, and is nondecreasing in its variables, except for the domain Q. It should be noted 
that the letter C{. . .) may depend on the physical parameters and the dimension N in some 
places, however we usually omit this dependence for simplicity. 

Our existence result of large weak solutions for fll.ip - fll.5p reads as follows. 
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Theorem 1.1. Let N = 2 or 3, Q G be a bounded domain of class C^'" with a G (0, 1), and 
the initial data po, uiq, do satisfy the following conditions: 

PoeL^iQ), Po > a.e. inQ, (1.8) 

moGL^+i(fi), mol{po=o} = a.e.inQ, l{po>o} G L (fi), (1.9) 

Po 

do G H2(fi), Idol < 1 in Q. (1.10) 

Then, there exists a constant eo := eo(A^, ^^) < 1 depending on N and Q (but independent of the 
physical parameters in U.l]) - fT7^) and the initial data), such that if do]\f := (ioAf(x) (the N-th 
component of do{x.)) satisfies 

l-doN<eo, (1.11) 

the initial- boundary value problem l^l.l\) - fl~S\) has a global weak solution (p, v, d) on I = It for 
any given T > 0, with the following properties: 

(1) Regularity: 

0<p a.e.mQr, p e C\lLl^^^{n) r\C\lL^m n L^^^Qt), 



V G L^(/,HJ(^])), pv G L-(/,L^(fi)) ^c°(/,L;:,^,(^])), 

\d\<l a.e. inQr, d G ^2(1, H2(fi)) n C°(J, Hi(r])), dtd & L^I {^1)) , 
(n • Vd)|gn = m the sense of trace for a.e. t G /, 

where p G [1,7), and r] G (0, (27 - A^)/A^). 

(2) Equations U.l\) and hold in {'D'{Qt))^~^^ , and equation U.3\) holds a.e. in Q^. 

(3) Equation U.l\) is satisfied in the sense of renormalized solutions, that is, p,v satisfy 

dth{p) + div[6(p)v] + [ph\p) - h{p)] divv = zn P'(M^ x J), (1.12) 
provided (p,"v) is prolonged to be zero on \ Vt, for any h satisfying 

b G C°[0, 00) n C\0, 00), \b'{s)\ < cs"^«, s G (0, 1], Ao < 1, 
and the growth conditions at infinity: 

\b'{s)\ < cs^\ s > 1, where c> 0, < 1 + Ai < 

(4) Regularity estimates: 



27 



:i.i3) 



sup(^(t) + ||d(t) -e^||L2(n)) + ||(Vv, VM)||l2(q,) + ||Vd||L4(Q^) < C{So,T,n), 
where Sq := S{0) = S{pq, mo, do), 62 = (0, 1), 63 = (0, 0, 1), and we have defined 

In particular, if Q is a ball Br := {x G M"^ | |x| < R} with R> 1, then the above constant 
eo can be chosen to be independent offl for any R> 1. Moreover, the constant C{Sq, T, Q) 
in fll.lSp can be replaced by a constant C(£^o, ^, ||do — eiv||L2{n)) independent offl. 



(5) In the case of N = 2, if, in addition, |do| = 1, then the weak solution satisfies |d| = 1, and 
the following finite and bounded energy inequalities: 

^ + T{t)<0 znV'il), 

(JjL 

S{t)+ [ J^{s)ds<So for a.e. tel. (1.14) 
Jo 

Remark 1.1. The proof of Theorem 11.11 remains basically unchanged if the motion of the fluid 
is driven by a bounded external force, i.e., when the momentum equations (11. 2p contain an 
additional term pf(x, t) with f being a bounded and measurable function. We remark here that 
we do not require any smallness condition on f . However, we are not clear whether the above 
theorem still holds with non-homogenous boundary condition in place of Neumann boundary 
condition "(n ■ Vd)!^^ = 0". In the proof of Theorem II. H we use the Neumann boundary 
condition only in order to deduce a maximum principle on d. 

Remark 1.2. We mention that the regularity requirement "do(x) G H^(ri)" is not optimal, for 
example, if we have "do(x) G W^'^(i7) with p > N" , then the above theorem still holds, and 
this can be shown by a standard approximate approach. On the other hand, we do not known 
whether "do(x) G H^(f2)" is the lowest regularity requirement, since it involves the problem of 
the Sobolev maps between two manifolds with the lower boundedness condition (II. lip . 

Remark 1.3. In view of the above regularity estimates in a ball, we can make use of a domain 
expansion technique to obtain a similar existence result of global weak solutions to the corre- 
sponding Cauchy problem, for which the expression of energy S{t) should be written in a form 
around some equilibrium state (poo, Vqo, gat) with poo > to make the energy integral sense (see 
0, Theorem 1.2]). Of course in this case, we can also establish a similar existence result of 
global weak solutions to the corresponding incompressible problem. 

We now describe the main idea of the proof of Theorem 11.11 For the Ginzburg-Landau 
approximation model to (ll.ll) - (ll.3p . based on some new estimates to deal with the direction field 
and its coupling/interaction with the fluid variables, Wang and Yu in S^] adopted a classical 



three-level approximation scheme which consists of the Faedo-Galerkin approximation, artificial 
viscosity, an artificial pressure and the celebrated weak continuity of the effective viscous flux to 
overcome the difficulty of possible large oscillations of the density, and established the existence of 
weak solutions. These techniques were develo ped in 22[] and p, ll2j for the compressible Navier- 



Stokes equations, we refer to the monograph [26| for more details. In the proof of Theorem 11.1 



we also adopt the three- level approximation scheme, so the key steps are to deduce the a priori 
estimates and to construct approximate solutions to the third approximate problem. Compared 
with the Ginzburg-Landau approximation model in [sol], however, the system (ll.ip - (ll.3p is much 
more difficult to deal with, due to the supercritical nonlinearity |Vdpd in (II. 3p . Consequently, 



not like that in [30|, one can not deduce the (sufficiently) strong estimate V d G L (J,L (fi)) 
directly from the basic energy inequality (I1.14p . Recently, Ding and Wen obtained the global 
existence and uniqueness of strong solutions to the 2D density- dependent incompressible model 
with small initial energy and positive initial density away from zero in j^, where they got V^d G 
L^(J, L^(f2)) from the basic energy inequality under the smallness condition of the initial energy. 
In fact, they first deduced Vd||L2(Q)/2 < So and iy6\\Ad\\l2(^Q^) < So + '^^ll I Vd| 1^4^^^^ by 
employing the basic energy inequality, and then made use of the inequality 

||Vd||t4(^) < C(fi)(||VM||^.(^)||Vd||^,(^) + ||VdrL.(f,)) 
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for some constant C{Q) depending on f2 C M^, which follows from the elliptic estimates and an 
interpolation inequality (see js], Lemma 2.4]), to infer that 

1|VM||^,(q^^<C(£o,||VMo1|l^(c)), (1.15) 

provided that the initial energy is sufficiently small. Motivated by this study, Jiang et al jl3] 
established the global existence of weak solutions to the corresponding compressible problem. In 
the current paper, we shall use another version of the interpolation inequality 

l|Vd||^4(f^) < C(l])(||VM||^2(f^)||d - ejvllL-(Q) + l|d - e7v||t4(Q)) 

for some constant C{Q) depending on i7 C {N = 2 or 3), from which the estimate fll.lSp can 
also be deduced for Q C if ||d — e7v||Loo(Q) is sufficiently small. Now the question is whether 
the smallness of ||d — eAr||Loo(Q) is guaranteed by smallness of the initial data ||do — GN\\i^oo(^Qy 
Fortunately, this is the case by applying the maximum principle to nonnegative lower bounds of 
solutions to the equations f 1 1.3 1) and the condition |d| < 1. Consequently, we deduce the desired 
energy estimates on d from the energy inequality. With these estimates in hand, we can adopt 
and modify the three-dimensional approximation scheme approach to show Theorem II. H if we 
can construct a solution to the following third approximate problem: 

dtp + dw{p\) = eAp, (1.16) 
9td + V ■ Vd = ^(Ad + /.(iVdHd), (1.17) 



/ (pv)(t) ■ *dx - / mo ■ *dx 

Jn Jn 



Jn 



/iAv + {p + A) Vdivv - VP - 6Vp'^ - e{Vp ■ Vv) (1.18) 

vdpr 



div(pv (g) v) - i/div Vd (g) Vd 



■ *dxds. 



where the tt,- dimensional Euclidean space X„ will be introduced in Section ^ e, S, (3 > are 
constants, and the smooth function felx) > satisfying 

r < fix) < 1, 

feix)=xiiN = 2; I f{x) = e-^iix>e-\ if iV = 3. (1.19) 

t < a; - fe{x) as £ 0, 

It should be noted that the third approximate problem above still enjoys the desired energy 
estimates (see Proposition 12.11) . thus it is easy to establish the unique solvability of the third 



approximate problem in the 2D case by following the same proof as in [10|. However, the proof 
in 



lo| can not be directly applied to the 3D case, and the difficulty lies in that we could not 
deduce a global estimate on ||9id||j;^oo(-7 l2(q)) (see 04.261) ) for the 3D approximate problem fl4.ll) - 
(14.61) . To overcome this difficulty, we introduce the cut-off function (11.191) to get a global estimate 
of II <9td|| 2,00(7 L2(Q)). On the other hand, we have to pay the price for this, namely, for the 3D 
approximate problem fll.l7p based on a cut-off function with Neumann boundary condition, we 
can not show |d| = 1 when |do| = 1. This is the reason why the solution in Theorem 11.11 does 
not satisfy |d| = 1 in three dimensions. 

The rest of paper is organized as follows. In Section [2] we deduce the basic energy equalities 
from the third approximate problem and derive more energy estimates on d under the assumption 
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fll.lip . In Section [3] we introduce the strong solvability of sub-systems in the third approximate 
problem, while the unique solvability of the third approximate problem is established in Section 
m Finally, we briefly sketch how to use the standard three-level approximation scheme to prove 
Theorem 11.11 in Section [5l 

2. A priori for the third approximate problem 

This section is devoted to formal derivation of the a 'priori energy estimates for the third 
approximate equations: 

+ div(pv) = eAp, (2.1) 
9t(pv) + div(pv ® v) + VP + 5Vp^ + e(Vp ■ Vv) 

= yuAv + (p + A)Vdivv - z/div(^Vd © Vd - i|VdpI^, (2.2) 

9id + V ■ Vd = 0(Ad + /,(| VdHd), (2.3) 

in a bounded domain C with initial data 

p(x, 0) = po > 0, d(x, 0) = do, v(x, 0) = vo, (2.4) 

and boundary conditions 

Vp-n|af, = 0, v|af, = 0, (n-Vd)|an = 0. (2.5) 

The a priori estimates will play a crucial role in the proof of existence. We consider a classical 
solution (p, V, d) of the initial-boundary problem fl2.1l) - fl2.5p with p > 0. 

2.1. Basic energy estimates 

We first deduce some basic energy estimates without any smallness condition imposed on the 
initial data. 

2.1.1. Maximum principle on \d\ 

The macroscopic average of the nematic liquid crystal orientation field d satisfies 

|d| < 1 in Qt, if |do| < 1 in C M^. (2.6) 

Next, we give a proof of (12. 6p for the reader's convenience. Multiplying the d-system (11. Sp by d, 
we obtain 

+ Iv ■ V|d|2 = ^^(Ad ■ d + /,(|Vd|^)|d|2). 

From the identity A|dp = 2|Vdp + 2 Ad • d it follows that 

dt{\d\^ - 1) - ^A(|d|2 - 1) < -V ■ V(|dp - 1) + 20|Vd|2(|d|2 - 1). (2.7) 

Now, letting d = \d\^ — 1 and d+ = maxjd, 0} > 0, multiplying ( 12. 7p by d^^ and integrating over 
Q, we integrate by parts and use the boundary conditions to infer that 

^ I d\d^< I (4e|Vdp + divv)c/^dx < ||4e| Vdp + |divv| ||Loo(f^) j d^dx. 
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Assuming that (v, d) satisfies the following regularity 

1146*1 Vdp + \diYv\\\L^i,L°°{n)) < oo, 

we are able to apply Gronwall's inequality to get fl2.6p immediately. 

We shall see that all the couples (v„, d„) in the third approximate solutions constructed in 
Section H] satisfy the regularity required above. We remark that (12. 7p becomes an equality in the 
two-dimensional case, and one can get by directly multiplying (12. 7p with |dp — 1 that 

|d| = lingr, if |do| = 1 in ^] C (2.8) 

2.1.2. Energy inequality 

Integrating by parts and utilizing the boundary conditions, one easily sees that the system 
(I2.ip - (l2.2p satisfies the energy conservation: 



^Esit)+ [ (/i|Vv|2 + (A + /x)|divv|2 + £<5/3p/5-2|Vp|2 + A£7p^-2|Vp|2)dx 
= -y I (Vd)'^Ad ■ vdx, 

where (Vd)^Ad := {didj) ATxivAd and 

^ , , [ fl Imp Ap< 5 b\ . 

A7V2^o~ ''^°^^7^^^^ / with m = pv. 



(2.9) 



Multiplying (12. 3p by —Ad and integrating by parts, we have 
i 

2 St 



/ |Vd|Mx + / lAdPdx 



/"(v Vd)- Addx + /" [/,(|Vd|2)|Vdp + 2 f'e{\^d\'){d,A,Vd,A,-V)d-d\d^ 

\ l<i,j<N ) 



< / (v Vd) ■ Addx + ^ / (|Vd|^ + Co|VdnvM|)dx, 
Jq Jn 

which, together with (12. 9p . implies 



-Ssit) + ||/i|Vv|2 + (A + p)|divv|2 + ue\Ad\^ + eS(3p^-^\Vp\^ + Ae-fp^-'^lVpl^Lnn) 



d 
Jt 

<9u [ (|Vd|^ + Co|VdnvM|)dx 



(2.10) 



where 



Ssit) := Ssip, m, d) := + ——p^ + -^pf^ + dx. 

For the two-dimensional case, recalling |d| = 1 and fe{x) = x for x > 0, we can deduce the 
standard energy equality. In fact, multiplying the equations (12. 3 p by Ad + |Vdpd, integrating 
by parts and using the boundary conditions, one deduces that 



1 d 

2di 



Vd|Mx + / |Ad+ |Vdpdpdx= /(vVd)-Addx 



n 



which, together with fl2.9p . yields 

/ [/x|Vv|2 + (A + /i)|divv|2 + ue{\Ad + |Vd|2d'2^ 
Jn 



whence 



+ eSPpf-^lVpl' + A£7p^-2|Vpndx + j^Ssit) = 0, 



t 

' [p|Vv|2 + (A + /i)|divv|2 + u9{\Ad + iVdpdH 
-'n 



(2.11) 



(2.12) 



+ e6l3p''-'\Vp\' + Ae^p^-'\Vp\']d^ds + £s{t) = £s,o, 
where £"5,0 := ^^^(Po, mo, do). 



2.1.3. Maximum principle on lower hounds 

The system (12. 3p possesses the following maximum principle on nonnegative lower bounds: 

di > dj^i if d^i > dQ^> for any given constant c?g-, (2.13) 

where 1 < i < N, and we have denoted the z-th component of d and do by rfj and doi, respectively. 
This conclusion will play a crucial role in this paper, so we give its proof here for the reader's 
convenience. 
Letting 

uJi = di- ^0-, u~ = min{a;i, 0} < 0, 
we can deduce from (12.31) that 

dtUi - eAuji = efe{\Vd\''){ui + dj - V ■ Voji. (2.14) 

Multiplying ( I2.14p by , and using the Neumann boundary condition, (I1.19P and Holder's 
inequality, we find that 

whin) + (^W'^^i \\h(n) 

= [ [^/,(|Vd|2)(u;, + rfo,.)-v Vu;.]u;rdx 
Jn 

= [ eM\Vd\'\)\ujr\M^ + l [ divv|u;r|2dx+ / eM\Vd\')d,^uj-d^, 
Jn ^ Jn Jn 

< {pl^^l'^Wi-^in) + ^l|divv||L-(n))||wrlli2(n) + ^^|Vd|2doi^-dx, 
which, together with the fact dQiUJ~ < 0, yields 

j^W^iWhin) < (2||^|Vd|2||^^^^^ + ||divv|U^(j,)) 
Hence, if we apply Gronwall's inequality to the above inequality, we obtain 



l-rWlli.(.) < ||a;-(0)||i.(,)e^«^^ll^l^^l^ll-(n)^"'^'-"-<"))'^^ = 0, 



which gives fl2.13p . 



2.2. More estimates under the small oscillation condition imposed on d 

To obtain more estimates on d under the small oscillation condition, we first introduce the 



well-known Nirenberg interpolation inequality (see 25|, Theorem]): 



Lemma 2.1. Let u belong to L'^(M^) and its derivatives of order m, V"^u, belong to L'"(M^), 
1 < r < oo. Then for the derivatives V^u, < j < m, the following inequality holds. 



^^u\\lp(^n) < Co\\V"^u\\lr(^N^\\u\\ "^f^y (2.15) 



where 



for all a in the interval 



- = - + a + 1 - a 



p n \r n J q 



— <a<l 



m 

(the constant Cq depends only on n, m, j, q, r, a), with the following exceptional cases: 

(1) If j = 0, rm < n and q = oo, then we make the additional assumption that either u tends 
to zero at infinity or u & L^(]R^) for some finite q > 0. 

(2) If 1 < r < oo, and m ~ j — n/r is a non-negative integer, then Ii2.15\) holds only for a 
satisfying j /m < a < 1. 

In addition, for a bounded domain Q (with smooth boundary) the above assertions hold if we 
add to the right side Ii2.15\) the term 

C{n)\\u\\L^n) 

for any q > I. All the relevant constants thus depend also on the domain. 

Next, we derive more estimates on d under the assumption that the initial value of liiv satisfies 

1 — Co < doN < 1, for some eo G (0, 1]. (2.16) 

It should be noted that the constant C (fl) in the following deduction will denote various positive 
constants depending on its variable fl, but the constants Co and Ci{Q)-Cs{Q) are fixed. 
First, one gets from the maximum principle that 

1 — Co < c/Ar(x, t) < 1 for any t > and any x 6 il. (2-17) 
Recalling |d| < 1, that is YliLi < 1, one finds that 



\di\ < ^1- d% < iorl<i<N-l, 
which combined with (12.17^ leads to 

||d - e7v||L-'(n) < Co^/e^. (2.18) 

Thanks to Lemma 12. H we have 

||VM|U4(f,) < C{n){\\V\\\l2f^^^\\u\\l^^^^ + ||m||l4(c)) for C with TV = 2 or 3, 
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which yields 



in))- 



(2.19) 



To bound the right hand of (12.191) . we shall use the following elliptic estimate: There exists 
a constant Ci(fi), such that 



VM||l2(^) < C'l(^])(||Ad||L2(^) + ||Vd||L2(f,)) for any Vd G H\n) 



(2.20) 



where Q C M^, which can be deduced from [26|, Lemma 4.27]. Thus, putting (I2.18I) - (I2.20I) 
together, we conclude that 



Vd||^4(f,) <C2{n)Coeo\\Ad\\l^^-^ + Cm\\d - e;v||io.(f,)||Vd||^ 



<C2in)Coeo\\Ad\ 



HQ) 
2 

L2(n) 



HQ) 



|d - eN\\i4(^a)] 



C(fi)(||Vd||^.(^) + ||d-e^||^.(^)), 



(2.21) 



where the constant C*2(fi) > 2 ^ only depends on Q. Utilizing (I2.2ip . ( l2.2Up . and Cauchy's and 
Holder's inequalities, we can deduce from fl2.10p that 



+ (A + /i)|divvp + u9\Ad\'' + e5(3p^'^\Vp\'' + Ae-fp^-''\V p\ 



|L2(Q), 



< C'3^i/(||Vd||^4(f,) + ||Vd||?.4,oJ|VM| 

< c^ev fc(^])C'oeo|| Ad||2,(^) + c(n)(||vd||^,(^) 



|d - eAr||L2(!^)) 



+ C'i(r])(y^(72(f^)C'oeol|Ad||L.(f,) + v/CXl^(||Vd||L2(f,) 
x(||Ad||L2(f7) + ||Vd||L2(f7)) 



|d - CatIIl: 



(2.22) 



CiiSl)^ C^mC + 2C2(f^)Co ) v^||Ad||^^.(^) 



1 



4C, 



l^d||L2(n) 



+ C(^])(||Vd||^2,f,, + ||d-e^|'2 



"T iiiJ - «Af|lL2(n)^ 

where C3 > 2~^ denotes a constant. Now, choosing eo := eo(^^) G (0, 1] such that 



(2.23) 



we get then 



/i|Vv|^ + (A + /i)|divv|^ + ylAdl" + ed/Spf-'lVpl' + Ae-fp^-'\Vp\ 
|d-e^,||?.,o^) < C(n)(£5(t) + 1), 



Li(n) (2.24) 



< C(n)(||Vd||L2(f^) - UNllL^in)) 

which, together with Gronwall's inequality, implies that 

Ssit) < C{Ss,o, T, Q) for any t e (0, T). 
Consequently, we can further infer that 



/i|Vv|' + (A + /i)|divv|' + ^|Ad|^ + e5l3p^-''\Vp\' + Ae^p^-''\Vp\ 



n7-2| 



(2.25) 



<C(^5,o,T,fi), 



11 



where Qt := (0,t) x Q for any t G (0,T). Moreover, from fEl^ . flT^ and (IT^ we get 

l|Vd||i4(,,L4(f,)) + ||VM||L2(n) + ||Vv||i.(,^L2(^)) < C{Ss,o,T,Q). (2.26) 

Finally, using (12 .26 p . Holder's and Sobolev's inequalities, we find from the equation (12.31) that 

ll^td||L4/3(j,L2(f^)) <C(||v ■ Vd||i4/3(j,L2(n)) + ll^d + A(|Vdnd||i4/3(j,L2(n))) 27) 
<C{Ss,o,T,Q). 

Similarly, we can also deduce that 

\\dtd\\L2ii,iHHn)r) <C{Ss,o,T,n), (2.28) 

where (H^(r2))* denotes the dual space of H^(fi). 

In addition, when fl = Br with R > 1, we can show that all the previous estimates on 
(p, V, d) are independent of Br, except for dtd. In fact, using (I2.19p and (I2.20p for VL = Bi, and 
scaling the spatial variables, we can obtain 



|Vd||^4(^^)<Co(||VM||^.. 



^n\\-Lo^{Br) 



^N\\-Li(^BR)) 



and 



|VM||l2(b,) < Co(||Ad||L2(B«) + ||Vd||L2(B^)). 



Hence, repeating the deduction process of (I2.24p . and employing the above two inequalities, one 
can have the following estimate: 



d 
dt 



Ss{t) + /i|Vv|^ + (A + /i)|divv|' + ^|Ad|^ + e6f3p('-^\Vp\'' + Ae-fp^-''\Vp\ 



o7-2| 



<^(l|Vd||^. 



(Q) 



|d - eAr||L2(n)) 



LHn) (2.29) 



where the constant C is independent of = Br for any R> 1. On the other hand, using (I2.18p . 
(I2.23P and Cauchy-Schwarz's inequality, we see from (12. Sp that 



1 d 
2rft' 



d-eiv||L2(iJH) + ^ll^d||2,,„_, =0 



(Br) 



Br 



/,(|Vd|2)(l-rf^) + l|d-ejv|Mivv 



dx 



I di vv I 



lHBr)\ 



GnWl^Br) 



< -||Vd||?2rR ^ , 
— 2 II ^^l,^(Br) 2 



-||divv||i2(B^) + C(0,p ^)\\d - eNWhf^BR)- 



Adding the above estimate to (12.241) . we get 



/i|Vv|2 + 0|Vd| 



+ (A + p)|divvp + ^|Ad|' + eS(3p^-^\Vp\'' + Ae7p^-^|Vp| 



Li(Q) 



+ ^ (Ssit) + ||d(t) - e^.||^,(^^)) < C (^Ssit) + ||d(t) - e^vl 



which, together with Gronwall's inequality, yields 



p|Vv|2 + 0|Vd| 
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+ (A + /i)|divv|^ + y |Ad|^ + e6Pp^-'\Vp\' + Ae^p^-'\Vp\ 

+ Ss{t) + \\d{t)-eN\\l2^BR,)<C{Ss,o,T,\\do-eN\\-L2^^^), for all tel. 
Summing up the above estimates, we conclude that 
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Proposition 2.1. Let N = 2 or 3, and Q C be a bounded domain and cq := eo(A^, ^) G (0, 1] 
satisfy (12.231) . Then the initial-boundary value problem enjoys the following a priori 

estimates, provided the initial data doN satisfies 1 — do^ < eo and |do| < 1. 

|d| < 1, in particular, |d| = 1 if |do| = 1 for N = 2, (2.30) 

||d — e7v||L°°(n) < C'ov^' (2-31) 

sup ||(d(t) - e^, Vd(t))||L2(f,) + ||V2d||L2(Q^) + ||Vd||L4(Q^) < Ci£s,o,T,n), (2.32) 
tei 

sup(||Vpv(t)||L2(^) + ||p(t)|U.(n)) + \\^^LHi,LHn)) < Ci£s,o,T,n), (2.33) 
tei 

l|f^td||L4/3(/,L2(n)) + ||c^td||L2(/_(Hi(f7)).) < C{£sfi,^,T). (2.34) 

In particular, if Q = Br with R > 1, then eo can be chosen to be independent of the domain 
Q = Br for any R> 1, and the constant C{Ss,o,T,Q) in (12.321) and (I2.33P can be replaced by a 
constant C{£sfi,T, ||do — eAr||L2(Q)) independent of Br. 

3. Strong solvability of sub-systems in the third approximate problem 

Before proving the unique solvability of the third approximate problem, we shall introduce 
two preliminary results. The first result is concerned with the global solvability of the Neumann 
problem on the equation (12.11) for given v. 

Proposition 3.1. Let < a < 1, Q be a bounded domain of class C^'", and po •= Po(x) satisfy 

po e W^^•°°(^]), 0<p<po<p<oo. 
Then, there exists a unique mapping 

y,,:L^ii,wl'^m^ C'{I,H\n)), 

such that 

(1) S^pQ ( v) belongs to the function class 

:= {p\p(iL^[l,W'^'''(Sl))^C''{l,W^^''{Sl)), dtpeL\l,L''iQ))}, 1< g < oo. 

(2) The function p = ^po(v) satisfies the following initial-boundary problem: 

dtp + div(pv) = eAp a.e. in Qt, 

p(x, 0) = Po a.e. in Q, (3.1) 

Vp ■ n\gn = 0, in the sense of traces a.e. in I. 

(3) ^poi^) is pointwise bounded, i.e., 

pg-Jo Mlwl,-(n)dr ^ ^^^^(v)(x,t) < pe^O l|v(r)|lwl,oo(o)dr^ ^ ^ ^^^^ ^ ^ ^ ^^3^2) 

(4) // ||v||icx>(j_wi>°°(n)) < i^v, then 

\\ypoi^)\\L^iU,mm < C(r])||po||^,i(f,)e^(-+'^")*, 

||V2^,„(v)||l2(qo < ^v^||po||Hi(n)«:.e^(-+«^)*, (3.3) 

||5t^Po(v)||L2(QO < C{n)Vi\\po\\HHn)K.e^'^'^+^'^^' 
for any t E I , where It '■= (0, t), and Qt := Q x It- 
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(5) ^poi^v) depends continuously on v, i.e., 

IK-^Pol^^l) ^ -^po('^2)](i)||L2(n) < G{Ky, e, T)t\\po\\ Hi {n)\\^l - V2||L°°(/t,Wi.°°(Q)), 

\\dt[^poi^i) - ^po(v2)](t)||L2(QO < G{K^,e,T)Vt\\po\\m{n)\Wi - vsllL-c/t.wL-CQ)) 

for any t G /, and for any \\vi\\L°°{i,wi-°°(ny) < Q'^'^ ||v2||loo(/ ■wi'°°(f7)) < ^^i;- The constant 
G is nondecreasing in the first variable and may depend on Q. 



Proof. Please refer to [26|, Proposition 7.39] and [10|, Proposition 3.1] □ 

The second result is on the local solvability of the Neumann problem for the system (12.31) 
with given v. 

Proposition 3.2. Let < a < 1, fl E be a bounded domain of class C^'", do := do(x) G 
H^in), K>0 and 



V G Yk 



:= h^iK, lldollH^m),^^) e (0,min{l,r}), 



lcO(/,H2(n)) ll'^t^llL2(/,Hi(n))^ 

Then there exist a finite time 

and a corresponding unique mapping 

^do "^Kiwith zn place of T) C\lj[ ,11^ {n)) , 
where hi is nonincreasing in its first two variables and Qt'^ := Q x := Q x [0,T^), such that 

(1) ( v) belongs to the following function class 

:= {d I d G L\lf, did G L\lf, (Hl(^]))*), 

(2) d = ^J^^(v) satisfies the following initial- boundary problem: 
(9fd + v- Vd = e(Ad + /s(|Vd|2)d) a.e. inQ^K, 

d(x, 0) = do a.e. in Q, (3.5) 

(n ■ Vd)|5f7 = in the sense of trace a.e. in , 

where fe>Ois defined by U.19\) . 

(3) ^d^(v) enjoys the following estimate: 

lK(v)||e. <C(ir,||do||H2(n),f^), 



where 



\f) l|2 J_ \\f) ||2 I II ||2 I II ||2 

1"^* llL2(/f ,Hi(n)) "t" ll"^* ,L2(n)) II llL°°(/Jf,H2(n)) II llL2(/f ,H3(f2)) 



Moreover (recalling Ii2.30\) and Ii2.13\) ). 

\d\ < 1, in particular, \d\ = 1 if |do| = 1 for N = 2; 
di{t) > rfg- for any t, if d^i > d^^ > where 1 < i < N. 
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(4) ^<^(v) continuously depends on v in the following s 



ense: 



I|[<(V1) - <(v2)](t)||L^(n) + ||V«(V0 - <(v2))||l.(q0 

< ViC{K, II do||H2(f7), fi)||vi — V2||loo(/^^l2(j^)) for any t E if- 

Proof. The above results on the two-dimensional non-homogeneous Derichlet problem have 
been shown under the higher regularity condition do E H^(f2) in 10|, where the appeared 
constants depend on || V^do||L2(n)- For the multi-dimensional Neumann problem considered here, 
we can obtain the same results by slightly modifying the proof in 0, Proposition 3.2] as follows. 

(1) We utilize a semi-discrete Galerkin method from 26|, Proposition 7.39] to construct a 
solution of the linearized Neumann problem, and a semi-discrete Galerkin method from [4] to 
construct a solution of the linearized non- homogeneous Derichlet problem in |10|]. 

(2) We make use of the elliptic estimate fl2.20p for the Neumann boundary condition to replace 
the corresponding elliptic estimate for the non- homogeneous boundary condition in 10| . This is 
the reason why the constant C in Proposition 13.21 is independent of || V^do||L2(Q)- 

(3) We use the following three-dimensional interpolation inequalities in Lemma 12.11 

< Cin){\\VMl2(^n)Ml2(^n) + Ml^u)) for u E H^{n), 

for u E H\n), 

(27l . Lemma 3.5]), 



\u\\L°°{n) 



\u\\LSn)<V2{\\Vu"'^ " 



\u\ 



(3.6) 
(3.7) 
(3.8) 



" \\u\\L^(n) 

\LHn)\\Hb{n) ioiuEH^in) (cf. 
to replace the following two-dimensional interpolation inequalities 

< C{n)\\u\\H2{n)\\u\\L2(n) for u E H^i^), (3.9) 
\H\h{n) <C{n)\\u\\Hi{n)\\u\\L2{n) for m e iJ^(fi), (3.10) 
MUn) < v^l|Vw||i2(^)||w||i2(^) for u E H',{Q). (3.11) 

The purpose of using fl3.9p - (l3.1ip above in l3| is to construct the term with a E (0, 1). When 
we use (I3.6p - (l3.8p in showing solvability of the three-dimensional problem, we can still construct 
the term with possible different a E (0, 1), with the help of Young's inequality in addition. 
Finally, we remark that the conclusions in Proposition 13.21 still hold in the 3D case, provided 

fe{x) =X. □ 



4. Unique solvability of the third approximate problem 

In order to obtain a weak solution of the initial-boundary problem f ll.ip - (ll.5l) . we first show 
the existence of solutions to the third approximate problem of the original three-dimensional 
problem ffTTl^-ffLHD: 



dtp + div(pv) = eAp, 

dtd + Y-Vd = e{Ad + /.(iVdHd), 



(pv)(t) • *dx - / mo ■ *dx 

Q Jn 

t 



Jn 



/iAv + (/i + A) Vdivv - AVp'^ 
- div(pv (g) v) - z/div ( Vd ® Vd - 



5Vp^ - eiVp ■ Vv 
VdPl 



*dxds 



(4.1) 
(4.2) 



(4.3) 
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for all t G / and any ^ G X„, with boundary conditions 

Vp-n\gn = 0, v|aQ = 0, (n-Vd)|an = 0, (4.4) 
and modified initial data 

pix,0) = poEW'^^in), 0<p<po<P<oo, (4.5) 

d(x, 0) = do G n'in), v(x, 0) = vo G X,, (4.6) 

where n denotes the outward normal to dQ, and e, 6, P, p,p > are constants. 



Here we briefly introduce the finite dimensional space X„. We know from [26 . Section 7.4.3] 
that there exist countable sets 

{Xi}Zi, < Ai < A2 < ■ ■ ■ , and 

C Wl^^iQ) n W2'P(n), 1 < p < 00, 

such that 

- /iA*i - {p + A)Vdiv*i = Ai*i, i = 1, 2, ■ • • , 

and 

I is an orthonormal basis in L^(f2) and an orthogonal basis in HQ(fi) with respect to 
the scalar product J^[pdjU ■ djV + (/i + A)divudivv]dx. We define a n-dimensional Euclidean 
space X„ with scalar product < ■, ■ > by 



Xj = span{^'j}"^;^, < u, V >= / u ■ vdx, u, v G X„, 

and denote by <^^„ the orthogonal projection of L^(f2) onto X„. 
4.1. Local existence 

With the help of Proposition 13.11 and 13. 2[ one can establish the local existence of a unique 
solution to the third approximate problem fl4.ip -( B^ by a contraction mapping argument. To 
this purpose, we rewrite the approximate momentum equations (14. 3 p as an operator form. 

Given 

p G C°(/,Li(fi)), dtpeL^inr), essinf p(x,t) >p>0, 

(x,t)gCT - 

we define, for all t E I, that 

X 

n 

by 

< ^p(t)^, w >:= / p(t)v ■ wdx, V, w G X„. 

Recall that all norms on X„ are equivalent, in particular, 

W^^'Pi(fi) and (Wo^'^2(r]))*-norms are equivalent on X„, (4.7) 

where (Wo^'^^(f]))* denotes the dual space of (W^^•P^(^])), ki and k2 are integers, and < k2 < 
00, 1 < < 00, < /ci < 1 and 1 < < 00 (or fci = 2, 1 < pi < 00). Note that this property 
of equivalent norms plays an important role in the estimates of velocity v. 
It is easy to see that 

^{x„,x„) < c{n) [ pit)dx, t G /. (4.8) 
Jq 
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On the other hand, we easily verify that -^^(4) exists for all t G / and 

||^,7,')IU(x„,x,o<^, (4.9) 

where <jSf (X„, X„) denotes the set of all continuous linear operators mapping X„ to X„. By 
virtue of (14.81) and f l4.9p . we have 

CiTl) — 

II-^pW-<=>iW^pwIU(x„,x„) < ^IIpiWIUhq), t e /. (4.10) 

Next, we shall look for T* C (0,Tf ] and 

V e A := {v G C(/:,X„) I G L2(/:,X„)}, /: := (0,r:) C (0,Ti^) 
with ||v||c7(j*,H2(n)) + \\dt^\\L2(i*,H^n)) < K for some K, such that 



/ ipy)it) ■ *dx - / mo ■ *dx 
Jn Jn 



Jn 



/iAv + (/i + A)Vdivv - AVp^ - 6Vp^ - e{Vp ■ Vv) (4.11) 

■ *dxds 



- div(pv (g) v) - z/div |^Vd (g) Vd 

for all t G [0, T„] and any ^ G X„, where p(t) = [^po{'v)](t) is the solution of the problem (13. ip 
constructed in Proposition 13. H d(t) = ^d^(v)(t) is the solution of the problem (13. 5 p constructed 
in Proposition [3]2l By the regularity of (p, d) in Propositions 13. II and 13. 2[ and the operator ^p(t), 
the equations (14. lip can be rephrased as 

v(t) = ((^mo + £ ^K(^,„(v), V, <(v))]d.) 

with mo = (pv)(0), where ^ := ^„ is the orthogonal projection of L^(i7) to X„, and 
^(p, V, d) = pAv + (p + A) Vdivv - AVp^ - ^Vp'^ - e{Vp ■ Vv) 

vdpr 



- div(pv (g) v) - z/div ( Vd (g Vd 

Moreover, one has 



dt^it) = ^^J^^^^^^^^^^^e4y^,i.)m [^^(^po(v),v,^do(v))](s)ds 

The authors in jiol . Section 4.3] have established the local existence of the problem (I4.ip - (l4.6p 
with non-homogenous boundary condition in place of the Neumann boundary condition for the 
system (14. 2 p by using a contraction mapping argument. In view of the proof in [13] and the 
previous preliminary results, we can immediately find that the results in [10|, Section 4.3] can be 
directly generalized to our problem (I4.ip - (l4.6p without essential changes in arguments. Thus, we 
have the following conclusion. 
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Proposition 4.1. There exist K := K{\\^mo\\x^, ||v(0)||x„, P~^) > and T* < T, such that 
T: A^A, r(w):=^[^(^)j|^mo + ^*[^^(^(w),w,^(w))](s)ds|, 

maps 

Bk,to = {w G a I ||w||c(I,^^x„) + \\dtW\\LHU„X^) <K}C C{lT,l,^n) 

into itself and is contractive for any < tq < T*, where one can take X„ = HQ(r2), T* has the 
form 

<T* = /i2(p, \\po\\m{n), WdoWn^n), P, K,T,n), (4.13) 
and /i2 is nonincreasing in its first three variables and nondecreasing in the fourth variable. 

Therefore, the map T possesses in Bk,t* a unique fixed point v which satisfies fl4.1ip . Thus, 
we have a solution (p = J^(v), v, ^(v)) which is defined in Qt* and satisfies the initial-boundary 
value problem fl4.1l) - fl4.6p for each given n. This means that we can find a unique maximal 
solution (p„,v„,d„) defined in [0,T„) x VL for each given n, where T„ < T. 

4-2. Global existence 

In order to show the maximal time Tn = T for any n, it suffices to derive uniform bounds for 
Pn, v„, d„ and ^n(PnV„). However, we need to impose an additional smallness condition as in 
(I2.16P to get the uniform boundedness of ||d„||L°°(/„,H2(n))- Fo^' simphcity, we denote 

(p, V, d, ^m) := (p„,, v„, d„, ^„(p„v„)). 

We mention that in the estimates that follow, the letters Gi{. . .), G2{. ■ •) and G{. . .) will denote 
various positive constants depending on its variables. 

First, one has the energy estimates as in Proposition 12. 1[ In fact, by virtue of the regularity 
of (p,v,d), we can deduce that (p, v, d) satisfies fl210|) for > 2, and fl241]) and fl212D for 
N = 2. Then, letting eo satisfy (12.231) and the initial data dg satisfy 

1 — eo{Q) < < 1 and |do| < 1, 

arguing in the same manner as in the derivation of fl2.30p - fl2.34p . we get 



|d| < 1, in particular, |d| = 1 if |do| = 1 for N = 2, 

||d - eAr||L°o(f7) < Go^/e^, (4.14) 

\\^/p^\\L^(Wn)) < G{Ss,o,T,n), (4.15) 

||v||L2(j„^Hi(n)) < G{£sfi,T,n), (4.16) 

||VM||l2(q,j + ||Vd||L4(Q^j + ||Vd|U^(,„,L2(f,)) < G{Ss,o,T,n), (4.17) 



||p||L-'(/„,L7(n) + l|f^td|lL4/3(/„,L2(n)) + \\dtd\\L2{i„,{n^n))") < G{Ssfl,T,Q). 

In particular, if f2 = with R > 1, the constant eo can be chosen to be independent of Q, 
and the above constant G{£sfi,T,Q) can be replaced by a constant G{£sfi,T, ||do — eiy\\i,2(^Q-j) 
independent of fl. 

With the help of fl4.15l) and fl4.16p . we can derive more uniform bounds on (p, v). Using 03.21) 
and fl4.16p . thanks to the norm equivalence on X„ (see (14. 7p ). we find that 

Gi{p,£s,o,n,T) < p < G2{p,Ss,o,n,T,Q), (4.18) 
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from which, f l4.15p and f l4.7p . it follows that 

l|v||co(/„,x„) < Gip,£s,o,n,T,n) (4.19) 

and 

ll^(pv)||co(j„,x„) < cH||pv||co(/„,L2(^)) < Gip,£s,o,n,T,n). (4.20) 
Applying KWh to ([33D, one gets 

\\p\\cOii„,m{n)) < G{p, \\po\\m{n),£5,o, n, T, fi). (4.21) 

Utilizing ( ICTl . (HT7D . ( HlSD - dOOll and (HZj), we obtain from gl2]) that 

||5tv||L2(j,^_x„) < G(p,p, ||po||Hi(n), ||do||H2(f7), ^5,0,^, T,fi). (4.22) 

Therefore, we have shown the uniform boundedness of || (p"\ p)||l-'(Qt„)' IIpIIco(J„,x„)' II v||co(/„,x,0' 
II (P v) 1 1 CO (/„ ,x„ ) and 1 1 9i v 1 1 ^,2 (/^^x„ ) • It remains to show the uniform boundedness of 1 1 d 1 1 ^oo (/^ (n) ) • 

Differentiating fl4.2p with respect to t, multiplying the resulting equations by dfd, recalling 
|d| < 1, we integrate by parts to infer that 

= 2 [ dtd- {OAdtd + 9fe{\Vd\^)dtd + 0/^(|Vdn(Vd : Vdtd)d - 9tv ■ Vd - v ■ Vdtd) dx 
Jn 

= ~29\\Vdtd\\l2^^) + 2 [ [(div9tv)d ■ dtd + dtv ■ Vdfd ■ d - v ■ Vdtd ■ dtd]d^ 

Jn 

+29 [ dtd ■ [f,{\Vd\')dtd + /^(|Vdn(Vd : V5id)d]dx 
Jn 

< -9\\Vdtd\\l^^^ + C(^^)(||9iv||^.(^) + ||v||^^(Q^)||9,d||^.(^) + ||9,d||^,(^^) 
+C{9){\\\dtd\'M\Vd\')U.^n) + |||a,dnVd|V;(|Vdr)|U.(o)), (4.23) 

where the last term on the right-hand side of fl4.23p can be bounded as follows. 

(1) The two-dimensional case: noting that fe{x) = x for the 2D case, we make use of Lemma 
12.11 and Holder's and Young's inequalities to see that 

C{9){\\\dtd\'M\Vd\')\\,.^n) + |||a,dnVd|Ve'(|Vdr)|U.(n)) 

< C(^)|||9,d||Vd|||i.(^) < Cmdtd\\l^n)\m\l^^^ 

< C{9,n){\\VdtdU.^n)\\dtdhHn) + l|5td||^,(^))||Vd||^,(^^ (4.24) 

< ^l|V9,d||^.(^) + ^^|^||9,d||^.(^)(l + ||Vd||^,(^)). 

(2) The three-dimensional case: recalling the definition of fe{x) in (ll.lQp . we use (13. 7p . and 
Holder's and Young's inequalities to get 

C{9){\\\dtd\'M\Vd\')\\,.^n) + |||a,dnVd|V;(|Vdp)|U.(^)) 



< 



c(^^)l|9,d||^.(^)(ll/,(|vdr)IU2(^) + |||vd|V;(|vdniU2(^)) 



< c(0,^])5-l(||v9,d|||,(^)||a,d||^,(^) + ||a,d||^,(^))||vd||L4(^) (4.25) 

< ^l|Va,d||^.(^) + ^(^'"^"''^) ||^.d||^.(^)(l + ||Vd||^.(^)) for eG(0,l). 
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Inserting and into (SSSD, we conclude that 

d 9 



-n — ^ iivdrL4(^) + iiv|i^.(Q^) + 1 + c{emv\\ 



which, by applying Gronwall's inequality, gives 

/ \ - 

\\dtd\\i.^^)<(\\dtd{miHn)+cm\dtA i2(/„,Hi(!^)) ) ^ 



-(ll^'l|lt4(Q^) + ll-llL-(Q^)+l) 



Noting that 

||5td(0)||L2(n) = ||^(Ado + |Vdo|Mo) - vo ■ Vdo||L2(n), 
we use firrgj) . (14:22|) and ( HlT]) to arrive at 

ll^td||ioc(/^_L2(Q)) + l|V^td||i2(/^_L2(Q)) <G{p,p, WpoWh^q), \\d.o\\n2{n), £5,o,n, e^'"^^ ,T,^l). 

(4.26) 

Recalling |d| < 1, from (14. 2p we get 

6^ [ |Ad|2dx<30^ [ |Vd|Mx + 3 [ |9td|Mx + 3 [ IvplVdpdx. (4.27) 
Jn Jn Jn Jn 

Similar to the derivation of f l2.2ip . the first term on the right-hand side of (I4.27P can be bounded 
as follows. 

36' [ |Vd|Mx < Se'C^mCoeoWAdWl^^^ + Cm\\^d\\l^n) + l|d - e^||^2(^)), (4.28) 
Jn 

where the constants 6*2(^2) and Cq are the same as in f l2.23p . Noting that C2{^)Co^o < i^Cs)^^ < 
4-1 by (1221, using f imp . fHTTp and f lOB]) . we find from (g^ZD and f lOg]) that 

^^|Ad|Mx <3^^C((])(||Vd||^.(^) + l) 

+ 3||5td||L2(o) + 3||v||^^(^)||Vd||22(^^ (^-29) 
<G{p,p, \\po\\m{n), \\do\\u'^n),£5fi,n,e^~'^ ,T,Q). 

Hence, by virtue of ^20^, KT7\f . ( KWf . flCT]) and the fact |d| < 1, 

||d||Lcx)(7„^H2(n)) <G{p,p, WpoWh^q), \\do\\m(n), £5,0, e^~'^^ fi). (4.30) 

The inequalities f Hl8|) . f Hl9|) . IK21\i . I K22\i and f lOOj) furnish the desired estimates which, 
in combination with Proposition 14. ![ give a possibility to repeat the above fixed point argument 
to conclude that T„ = T, and moreover, the global solution (p„, v„, d„) is unique. To end this 
section, we summarize our previous results on the global existence and uniqueness of a solution 
(Pn,v„,d„) to the third approximate problem (14. ip - (14. 60 as follows. 

Proposition 4.2. Let the constant eg > (depending on Q) satisfy (I2.23p . 

5 > 0, l3 > 0, e > 0, and < p < p < 00. 
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Assume that Q C is a bounded C^'"-domain (a G (0,1)), and the initial data (po,mo,do) 
satisfies 

l-doN<eo, |do| = l, doGH2(l]), (4.31) 
0<p<Po<P, PoeW''°°{n), voGX„. (4.32) 

Then the third approximate problem ( [y^.i[ j-( [775| ) possesses a unique triple {pn, v„, d„) with the 
following properties: 

(1) Regularity: pn satisfies the same regularity as in Proposition \3.1[ v„ G C°(/, X„), (9tV„ G 
L^(/, X„), d„ satisfies the same regularity as in Proposition \3.2\ with T in place of . 

(3) (p„,v„,d„,) solves (JjJ^ o^^'^d o,-^- Qt, and satisfies ( [^.3| ) and (pn, v„, d„)|t=o = 

(po, vo,do). 

(3) Finite and bounded energy inequalities hold in the 2D case: 

j^82{t) + ^"(t) + j^e5l3pi~-^\Vpn\\t)di^ < m V'{I), (4.33) 

and 

£sit) + £ (^"(s) + ^£5/3p^'|Vp„|2(s)dx^ ds < ^,(po,mo,do) := £s,o (4.34) 

a.e. in I, where J'"(t) := J'(p„, v„, d„) anc? £:^(t) := £^5(p„, m„, d„) with m„ = p„v„. 
f/^^ Additional uniform estimates: 

|d| < 1 in Qx, in particular, |d| = 1 if |do| = 1 /or N = 2, (4.35) 
||d„ - eN\\L^{n) < Co^/e^, (4.36) 
sup ||(d(t) - e^, Vd(t))||L2(n) + ||VM||l2(q^) + ||Vd||L4(Q^) < ^(^5,0,^), (4.37) 

tG/ 

l|f^td||L4/3(7-^L2(f^)) + ||9(d||i2(j^(Hi(n))*) < ^(£^5,0,^), (4.38) 
sup(||v/;^v(t)||L2(n) + ||p(t)IU7(n)) + ||Vv|U2(,,L2(f,)) < G(£5,o,^^), (4.39) 

y^||Vp„||L2(Q,) < G(£5,o,5,fi), (4.40) 
\\PnL^,^ <G{£s,o,e,5,n), (4.41) 

(see ^2^, Section 7.7.5.2] for the proof of ^4-4(^ (^nd ^4-4^^ ) j where G is a positive constant 
which is independent of n and nondecreasing in its arguments, and may depend on T . 
Moreover, if e is not explicitly written in the argument of G, then G is independent of e as 
well. 

(5) In particular, ifQ = with R> 1, then eo can be chosen to be independent of the domain 
Q = Bji for any R > 1, and the constant G in f l4.37p and f l4.39p can be replaced by a 
constant G{Ss,o,T, ||do — ejv||L2(Q)) independent of B^. 



5. Proof of Theorem 11.11 

Once we have established Proposition 14.21 we can obtain Theorem 11.11 by employing the 
standard three-level approximation scheme and the method of weak convergence in a manner 
similar to that in js], |22| for the compressible Naiver-Stokes equations. These arguments have 
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also been successfully used to establish the existence of weak solutions to other models from fluid 



dynamics, see the 2D problem of fll-ip -l lir^ in 10|], and the 3D Ginzburg-Landau approximation 



model to fll.ll) - fll.3p in [30| for example 



First we can construct a solution sequence v„,d„) by Proposition \A.2\ using the related 
uniform estimates in Proposition 14.21 and standard compactness arguments, we can obtain the 
weak limit (pj,Ve,de) of the solution sequence (p„,Vn,d„) as n — )■ oo, taking subsequences if 
necessary, which is a weak solution of the following second approximate problem: 

dtpe + div(p^vj - eAp^ = in T>'{Qt), 

dtipeVe) + djipey^evi) - pAv, - {p + A)VdivV, + V Ap] + 6V p^, 

f iVd Pl\ 

+ i/div Vd, ® Vd, - ' + eVp, ■ Vv, = in {V\Qt)Y , 



2 

dA, + V, ■ Vd, = ^(Ad, + /,(|Vd,nd,), a.e. in Qt. 

with boundary conditions 

Vpe -1113^ = 0, Velan = 0, (n ■ Vds)|an = 0, (5.1) 

and modified initial data f l4.3ip -( P:.32p . where 5 > 0, /3 > max{7,8}, and e > 0. Moreover, the 
solution (pe, V£,de) enjoys the finite and bounded energy inequalities (14. 330 - 04. 340 . and uniform 
estimates (l4:35D - fOTD . 

We proceed to utihze the related uniform estimates and standard compactness arguments 
to obtain the weak limit (p5,V5,d5) of the weak solution sequence (pe,Ve,de) to the second 
approximate problem as e — )■ 0, taking subsequences if necessary, which is a weak solution of the 
following first approximate problem: 

dtpf, + div(p5V5) = in V'{Qt), 

dtips^s) + dj{pswsvl) - fiAvs - (p + A)Vdivv5 + VAp] + SVp^s 
+ z/div (vds ® Vds - rZ^^ = in iV'iQT)f, 

dtds + vs ■ \/ds = e{Ads + iVd^^d^), a.e. in Qt 

with boundary conditions (15. ip and modified initial data (I4.3ip - (l4.32p . Moreover, the solution 
{ps,^^s,ds) enjoys the estimates (14.350 - 04.390 and inequalities (I4.33p - (l4.34p with e = 0. Here 
we remark that it is easy to verify the convergence of /e(|Vdep) to jVd^p as e — )■ in three 
dimensions. 

Using the uniform bounds given in (I4.36p -( l4.38p with d^ in place of d„, applying the Arzela- 
Ascoli theorem and Aubin-Lions lemma, and taking subsequences if necessary, we deduce that 

d, ^d5 Strongly in C'^{I,L\Q)) n LP{I,H\Q)) n L^{I,W^-''{Q)) (5.2) 

for any p G [l,oo), q G [1,6) and r G [1,2), which, recalling the definition of fs, implies that 
(taking subsequences if necessary) 

/sdVd^Hd^ ^ iVd^pd^ as£^0 for a.e. x G (5.3) 

Thus, using Vitali's convergence theorem, and recalling the uniform in e boundedness of || Vd£||L2(Qj,), 
we infer that 

/,(| Vd.Hd, ^ \Vds\^ds strongly in L^Qt) for any r G [1, 2), 

and weakly in L^((5t)- 
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In addition, we have the regularity G L^(J,H^(f2)) and dtds G L'^{I, (H^(f2))*). In view of 



26l . Proposition 7.31], we get consequently 

d5 G C°(/,Hi(fi)). (5.5) 

Finally, we can also obtain a weak solution (p, v, d) of the original problem fll.ip -f lTTil) with 
boundary conditions "v|aQ = and (n ■ Vd)|af7 = 0", and modified initial data fl4.3ip - p.32p . 
which is the weak limit as 5 — )■ of the weak solution sequence {ps, v^, d^) of the second approxi- 
mate problem. It should be noted that the modified initial energy in fl4.34p - fl4.39p can be further 
chosen to be independent of 6, in other words, the term £s,o in fl4.34p - fl4.39p can be replace by a 
positive constant Sq := supQ^g^i{£s,o}- Hence, the weak solution (p, v, d) enjoys the same esti- 
mates as in Theorem II. 1[ Applying an approximation argument to the initial data, the modified 
initial data fl4.3ip and fl4.32p can be relaxed to fll.SD - fll.lip . Consequently, we can obtain the 
desired Theorem II. 1[ We refer to [lol. Mi or 0, l26l | for the omitted details of the proof of the 



limit process and the renormalized solutions fll.l2p . 
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